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ABSTRACT 

Fully packed loop models on the square and the honeycomb lattice constitute 
new classes of critical behaviour, distinct from those of the low-temperature 0(n) 
model. A simple symmetry argument suggests that such compact phases are only 
possible when the underlying lattice is bipartite. Motivated by the hope of identi- 
fying further compact universality classes we therefore study the fully packed loop 
model on the square-octagon lattice. Surprisingly, this model is only critical for 
loop weights n < 1.88, and its scaling limit coincides with the dense phase of the 
O(n) model. For n = 2 it is exactly equivalent to the selfdual 9-state Potts model. 
These analytical predictions are confirmed by numerical transfer matrix results. 
Our conclusions extend to a large class of bipartite decorated lattices. 

PACS numbers: 05.50. +q, 11.25.Hf, 64.60.Ak, 64.60.Fr 
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1 Introduction 



Compact polymers, the continuum limit of random walks that are constrained to visit 
every site of some lattice £, are intriguing in so far as their critical exponents depend 
explicitly on £. Whilst first observed numerically |L|, this curious lack of universality 
was firmly established through the exact solution of the compact polymer problem on 
the honeycomb @, §] and, very recently, the square lattice j|, |J. 

However, not every lattice can support a compact polymer phase. To see this, consider 
more generally an 0(n)-type loop model defined on £, in which each closed loop is 
weighed by n, and each vertex not visited by a loop carries a factor of t. It is well- 
known that for |n| < 2 this model possesses a branch of low-temperature (t being the 
temperature) attractive critical fixed points || f|j with critical exponents that do not 
depend on £, even when £ is not a regular lattice but an arbitrary network ||. On 
the other hand, whenever the model is invariant under t — > —t, as is the case if £ 
can only accommodate loops of even length, this symmetry allows for a distinct zero- 
temperature branch of repulsive fixed points 0, with the n — ► limit representing the 
compact polymer problem. That the critical behaviour of this class of fully-packed loop 
(FPL) models depends on £ is readily seen from the solutions of the honeycomb and 
the square case given in Refs. || [4], Namely, the continuum limit of these models 
can be described by a conformal field theory (CFT) for a fluctuating interface, where 
the fully-packing constraint forces the height variable to be a vector, with a number of 
components that depends on the coordination number of the lattice at hand. 

This t — > —t symmetry argument, originally put forward by Blote and Nienhuis [|TJ, 
prompts us to conjecture its inverse: Whenever £ allows for loops of odd length, so 
that the symmetry is destroyed, the renormalisation group flow can be expected to take 
us to non-zero t, eventually terminating in the dense, universal O(n) phase. Support 
for this conjecture so far comes from numerics in the case of the triangular lattice |§, 
and recently for a class of decorated lattices interpolating between the square and the 
triangular lattices fTT)|j.p| 

2 Although belonging to the universality class of the square lattice FPL model {§ § the FPL model 
on the square-diagonal lattice does not constitute a very good counterexample, since the fully-packing 
constraint actually prevents the loops from occupying the diagonal edges. (Note that the proof given in 
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Accepting for the moment the validity of this conjecture however leaves us with an 
infinite set of bipartite lattices, each one being a potential candidate for a novel univer- 
sality class of compact polymers. This perspective is especially appealing in the light of 
the constructive point of view taken in Refs. [0, |, [|, P}| . In these papers new CFTs 
were explicitly constructed, based on purely geometrical considerations applied to the 
FPL model in question. On the other hand, if the bipartite lattices generate an entire 
family of distrinct CFTs, this gives rise to important classification issues. In particu- 
lar, the applicability of compact polymer models to the protein folding problem [12], [L3| 
implies that one would like to understand on which microscopic parametres (bending 
angles, coordination number, steric constraints) the resulting conformational exponents 
do depend. 

In this paper we examine FPL models on a class of bipartite lattices, in which every 
vertex of a regular (square or honeycomb) lattice has been decorated. A renormalisation 
group (RG) argument, essentially amounting to a summation over the decoration, reveals 



that the Liouville field theory construction mM should really be based on the undecorated 



lattice, but with bare vertex weights that depend on the loop fugacity n. This leads to a 
novel scenario in which, depending on n, the model may either renormalise towards the 
dense phase of the O(n) model or flow off to a non-critical phase, even for n < 2! 

The case of the square-octagon lattice, shown in Fig. [I], is investigated in detail. This 

lattice can be thought of as a square lattice in which each vertex has been decorated with 

a tilted square. Our interest in the square-octagon lattice stems from the fact that it is 

bipartite and has the same coordination number as the honeycomb lattice, but enjoys 

the symmetry of the square lattice. In particular it will enable us to assess whether the 

critical behaviour of compact polymers on a lattice £ depends only on its coordination 

number^, only on the bond angles, or on a combination of both these parameters. Our 

analysis suggests that the corresponding FPL model belongs to the dense O(n) phase 

for n < 1.88, whilst for n > 1.88 a finite correlation length is generated. For n = 2 we 

show rigorously that the model is equivalent to the (non-critical) 9-state Potts model. 

Ref. |0| is also valid for n^Q.) 

3 In the protein folding language this determines the number of close contacts per monomer of the 

folded chain. 
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Figure 1: Fully packed loops on the square-octagon lattice. In the corresponding transfer 
matrix, periodic boundary conditions are imposed across a strip of width L loop segments 
(here L — 4). The state space is that of all well- nested, pairwise connections amongst 
the L dangling ends in the upper row. 

The analytical results are confirmed by numerical transfer matrix calculations on strips 
of width up to L max = 18 loop segments. 

Having introduced the models in Section [| we present the analytical results in Sec- 
tion |3] and the numerics in Section Our results are discussed in Section |5|. 



2 The models 

A fully packed loop (FPL) model on a lattice £ is defined by the partition function 

^fpl = £ (2- 1 ) 

Gfpl 

where the sum runs over all configurations (?fpl of closed loops drawn along the edges 
of £ so that every vertex is visited by a loop. Within a given configuration a weight n is 
given to each of its N loops. 

An FPL model on £ can be generalised to an O(n) model by lifting the fully packing 
constraint and further weighing each empty vertex by a factor of t. Physically t corre- 
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Figure 2: Phase diagram of the 0(n) model on the honeycomb lattice. 

sponds to a temperature, the FPL model thus being the zero-temperature limit of the 
0(n) model. When £ is the honeycomb lattice, the resulting phase diagram is as shown 
on Fig. Q [|IJ. For \n\ < 2 three branches, or phases, of critical behaviour exist. Since £ 
is bipartite, the resulting t — > —t symmetry allows for a compact phase at t = 0, [3J, 
as discussed at length in the Introduction. For t > Nienhuis has found the exact 
parametrisation of a dense and a dilute phase, and determined the critical exponents as 
a function of n || . 

For our discussion of the square-octagon FPL model we shall need the corresponding 
parametrisation for the 0(n) model on the square lattice. The definition of the partition 
function is now slightly more complicated, since each vertex can be visited by the loops 
in several ways that are unrelated by rotational symmetry. An appropriate choice is 

Z Q{n) = J2 t Nt u N "v N *w N -n N , (2.2) 
G 

where N t , N u , N v and N w are the number of vertices visited by respectively zero, one 
turning, one straight, and two mutually avoiding loop segments. It is convenient to 
redefine the units of temperature so that t — 1. 

Nienhuis |7|, [L5| has identified five branches of critical behaviour for the model (2.2) . 
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The first four are parametrised by 



w. 



u. 



1 - 2 sin 



1 + 2 sin 



4^ c sin m cos [\--A, 



w. 



1 + 2 sin 



'0' 



?? 



-2cos(2#), 



(2.3) 



where 9 G [(2 — b)iv/2, (3 — b)iv/2] corresponds to branch b = 1,2,3,4. It has recently 
been noticed that the edges not covered by the original ('black') loops form a second 
species of closed ('grey') loops, each one occuring with unit weight ||TT|. Lifting the 
fully-packing constraint implies that the two loop flavours decouple, and each of them 
can independently reside in either of the two critical phases (dense or dilute) discussed 
above. The black (resp. grey) loops are dense on branches 2 and 4 (resp. 1 and 2), and 
dilute on branches 1 and 3 (resp. 3 and 4). On branches 1 and 2 the grey loops contribute 
neighter to the central charge, nor to the geometrical (string) scaling dimensions, and in 
the scaling limit these two branches are thus completely analogous to the dilute and the 



dense branches of the O(n) model on the honeycomb lattice fll 
The last critical branch, known as branch 0, has weights 



u. 



Wc 



0, -3 < n < 1, 



(2.4) 



and can be exactly mapped onto the dense phase of the 0(n+l) model J7|, or equivalently 
to the selfdual in + l) 2 -state Potts model 111. 



3 Renormalisation group analysis and an exact map- 
ping 

At first sight it would seem that the continuum limit of the FPL model ( |2.1| ) on the square- 
octagon lattice should be described by a Liouville field theory for a two-dimensional height 
field, since the lattice has the same coordination number as the honeycomb lattice ||. 
However, we shall presently see that only one height component survives when applying 

5 



n 



Figure 3: First step in a real-space renormalisation of the square-octagon lattice FPL 
model. The renormalised vertices get weighted as shown. 

the appropriate coarse graining procedure to the two-dimensional microscopic heights 
defined on the lattice plaquettes. 

Consider performing the first step of a real-space renormalisation group (RG) trans- 
formation of Eq. ( |2.1| ) , by summing over the degrees of freedom residing at the decorating 
squares. In this way the decorated vertices transform into weighted undecorated vertices, 
as shown on Fig. || The renormalised model is then simply the O(n) model on the square 
lattice ( |2.2| ), but with some particular 'bare' values of the vertex weights. Defining again 
the empty vertex to have unit weight, these bare weights read 



1 1 

u — — , v = 0, w = —. 
n n 



(3-12 



Following the standard procedure |16| microscopic heights can be defined on the lattice 
plaquettes by orienting the loops and assigning a vector, A, B or C, to each of the three 
possible bond states: A (B) if the bond is covered by a loop directed towards (away 
from) a site of the even sublattice, and C if the bond is empty. When encircling an 
even (odd) site in the (counter) clockwise direction the microscopic height increases by 
the corresponding vector whenever a bond is crossed. As was first pointed out in Ref. || 
the fully-packing constraint leads to the condition A + B + C = 0, whence the height 
must a priori be two-dimensional. However, the RG transformation that we have just 
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Figure 4: RG flow in the square-octagon FPL model. After tracing over the decoration, 
the bare value of 1/u is given by the dashed line. For n < 1.88 the flow is directed 
towards the attractive branch of dense fixed points, whilst for 1.88 < n < 2 the system 
renormalises towards the high-temperature disordered phase. 

applied lifts the fully-packing constraint, due to the appearance of the bottom left vertex 
of Fig. |^. Defining now the sublattices with respect to the renormalised (square) lattice 
we have the additional constraint 4C = 0, whence the coarse grained height field should 
really be one-dimensionalQ, and 0(n)-like behaviour is to be expected. Also note that it 
clearly suffices to define the microscopic heights on the octagonal plaquettes in order to 
obtain a continuous height field defined everywhere in IR 2 by the usual coarse graining 
procedure [0 . 



The reason that the renormalised FPL model is still interesting is that the bare vertex 



weights (|3.1| ) are now some fixed functions of the loop fugacity n, rather than arbitrary 
parameters that can be tuned to their critical values. This constitutes an interesting 
situation which has not been encountered before. We shall soon see that it implies that 
the FPL model (|2.1| ), unlike any other loop model studied this far, is only critical within 
a part of the interval |n| < 2. 

In Fig. |] we show l/u c , the weight of the empty vertex relative to that of a turning 



'See Ref. |Tl| for similar examples of such a reduction of the dimensionality of the height field. 
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loop segment, as a function of n for the critical branches 1 (dilute phase) and 2 (dense 
phase) of the O(n) model on the square lattice; cfr. Eq. Q2.3| ). In analogy with the 
honeycomb case the dense and dilute branches again consist of respectively attractive 
and repulsive fixed points. With the bare value 1/u given by Eq. (|3.1|) the subsequent 
RG flow must therefore be as schematically indicated on the figure. For n ~ 1.88 there 
is an intersection between the bare value and that of the dilute branch, and for n > 1.88 
we can therefore expect the flow to be directed towards the high-temperature disordered 
phase of the 0(n) model. In other words, a finite correlation length (roughly the size 
of the largest loop in a typical configuration) is generated and the model is no longer 
critical. 

Of course we should be a little more careful, since u is not the only parameter in the 
model. Whenever the bare weights Q3.1|) do not intersect one of the five branches of fixed 
points, v and w will flow as well. In particular, v will in general flow towards non-zero 
values, since the turning loop segments always occur with finite weight, and these are 
clearly capable of generating straight loop segments on larger length scales. The essential 
point is that for 1.88 < n < 2 empty vertices will begin to proliferate, and there is no 
physical mechanism for halting the flow towards the disordered phase. f] 

The point n = 2 merits special attention. Here the bare weights are 

u = w= 1 -, v = 0, (3.2) 



which coincides with the fixed point values on branch 0; see Eq. ( |2.4| ). Invoking Nienhuis' 
mapping [0] the n = 2 FPL model is therefore exactly equivalent to the selfdual 9-state 
Potts model, which is of course again non-critical |L7| . 



4 Transfer matrix results 

In order to confirm the analytical predictions given in Section |3| we have numerically 
calculated effective values of the central charge c and the thermal scaling dimension xt 
on strips of width L = 4, 6, . . . , 18 loop segments. To this end we adapted the connectivity 
basis transfer matrices described in Refs. J?], |] to the square-octagon lattice. The working 
5 The flow cannot be towards branch since this is a repulsive fixed point. 
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n 


c(4,8) 


c(6,10) 


c(8,12) 


c(10,14) 


c(12,16) 


c(14, 18) 


0(n) 


0.0 


-1.9784 


-1.9862 


-1.9880 


-1.9924 


-1.9963 


-1.9980 


-2.0000 


0.5 


-0.8898 


-0.8975 


-0.8729 


-0.8488 


-0.8338 


-0.8259 


-0.8197 


1.0 


0.0706 


-0.0856 


-0.0965 


-0.0701 


-0.0434 


-0.0246 


0.0000 


1.5 


0.9390 


0.6602 


0.5504 


0.5144 


0.5117 


0.5227 


0.5876 


2.0 


1.6484 


1.4844 


1.4495 


1.4362 


1.4225 


1.4068 


1.0000 



Table 1: Three-point estimates for the central charge, compared with exact results for 
the dense phase of the O(n) model. 

principle of these transfer matrices is illustrated in Fig. [TJ: To determine the number of 
loop closures induced by the addition of a new row of vertices it suffices to know the 
pairwise connections amongst the L dangling ends of the top row. For L even, the 
number of such connections is |7]] 

L/2 



(4.1) 



where c r , 



(2m)! 



. . are the Catalan numbers. Thus, the transfer matrix for a strip of 

m!(m+l)! ' 1 

width L has dimensions x a^, and a sparse matrix decomposition can be made by 
adding one site of the lattice at a time, rather than an entire row. The size of the largest 
matrix employed is given by ai& = 6, 536, 382. 

The effective central charge c(L, L + 4) has been estimated by three-point fits of the 
form 0, 

irr A 

(4.2) 



/ (L) = / (oo)-^ + 4 + 



QL 2 L 4 

applied to the free energy per site fo{L') with V = L, L + 2, L + 4. Similarly, effective 
values x t (L,L + 2) of the thermal scaling dimension were found from two-point fits of 
the form pOl § 

9.TTT, R 

(4.3) 



/ 1 (L)-/o(L) = y + - + - 



where fo(L) and fi(L) are related to the ground state and the first excited state of the 
transfer matrix spectra in the usual way. 

The numerical results are given in Tables [l] and 0. For n < 1.5 we see the expected 
convergence towards the exact values of the O(n) model in the dense phase, which read 
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n 


z t (4,6) 


z t (6,8) 


x t (8,10) 


x t (10,12) 


z t (12,14) 


^(14,16) 


x t (16,18) 


0(n) 


0.0 


1.2942 


1.9142 


2.5336 


2.5970 


1.9830 


1.9849 


1.9890 


2.0000 


0.5 


1.3792 


1.8833 


2.4173 


1.5522 


1.5201 


1.5557 


1.5685 


1.5843 


1.0 


1.2480 


1.2843 


1.2895 


1.2856 


1.2799 


1.2745 


1.2701 


1.2500 


1.5 


0.8593 


0.7842 


0.8088 


0.8392 


0.8664 


0.8876 


0.9030 


0.9482 


2.0 


0.5488 


0.5708 


0.5629 


0.4188 


0.3864 


0.3595 


0.3369 


0.5000 



Table 2: Two-point estimates for the thermal scaling dimension, juxtaposed with exact 
values for the dense O(n) model. 



El 



6e 2 



2e + 1 



1 - - 1 — > xt = , x (4-4) 
l-e 2(1 -e) V ; 

with e = - arccos(n/2). For n = 1.5 the convergence is rather slow, especially in the case 
of c, reflecting a large crossover length. 

As predicted by theory, the FPL model is no longer critical at n — 2. This is par- 
ticulary visible from the monotonic decrease of the Xt estimates, which are well below 
the exact O(n) value Xt = 1/2. For a system with a finite correlation length, £ < oo, 
the effective values for c should eventually tend to zero. The fact that we observe rather 
large effective values is in agreement with Ref. J7J, and rather predictable since £ is much 
greater than the largest strip width used in the simulations [p2 |. For comparison we 
performed similar computations for the 9-state Potts model in its loop representation 
finding again effective values of c in the range 1.3-1.4. 



5 Discussion 

Having seen that two of the simplest two-dimensional lattices (square and honeycomb) 
give rise to distinct compact universality classes, it would be tempting to conjecture that 
an FPL model defined on any new lattice leads to different critical exponents and has 
a new CFT describing its continuum limit. In the present paper we have demonstrated 
that this is far from being the case. Even within the very restricted class of bipartite 
lattices, fulfilling the t — > — t symmetry requirement, any lattice that can be viewed as a 
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decorated square or honeycomb lattice is likely to flow away from the compact phase by 
virtue of an RG transformation analogous to the one presented in Section [3[ 

Despite the curious lattice dependence of the compact phases, it thus appears that 
the number of distinct universality classes is very restricted. We recall that the con- 
tinuum limit of all loop models solved to this date can be constructed by perturbing a 
SU(iV)fc = i Wess-Zumino-Witten model by exactly marginal operators and introducing an 
appropriate background charge [Q. It would be most interesting to pursue the physical 
reason why only the cases N = 2 (the O(n) (|TTJ, Potts [^] and six- vertex models), 
N = 3 (the FPL model on the honeycomb lattice ||), and N = A (the two-flavoured 
FPL model on the square lattice [|, [5]]) seem to occur in practice. 

The square-octagon lattice FPL model studied here turned out to be interesting in 
several respects. First, it provides us with the first example of an non-oriented [p5| , 
bipartite 0, [1(J lattice for which the scaling properties of compact and dense polymers 
are identical. In particular, the exact value of the conformational exponent 7 is 19/16 



2% U, indicating a rather strong entropic repulsion between the chain ends. Second, 
the square-octagon model presents a novel scenario in which the same fully packed loop 
model may renormalise towards different conformal field theories, or even flow off to a 
non-critical regime, depending on the value of the loop fugacity |n| < 2. In particular one 
might be able to 'design' a decorated lattice with bare vertex weights that simultaneously 
intersect those of the dilute 0(n) phase for some value of n. This could be a starting 
point for gaining a microscopic, geometrical understanding of the Coulomb gas charge 



asymmetry |21| which was shown in Ref. |TTJ to distinguish between the dense and dilute 
phase of the O(n) model. Finally, our model proves that the scaling properties of compact 
polymers do not depend exclusively on either bond angles or coordination number, but 
rather on a combination of these two parameters. 
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